Deconfined criticality flow in the Heisenberg model with ring-exchange interactions 
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Quantum transition points in the J-Q model - the testbed of the Deconfined Criticality theory 
- and the SU(2)-symmetric discrete NCCP 1 representation of the deconfined critical action are 
directly compared by the flowgram method. We find that the flows of two systems coincide in 
a broad region of linear system sizes (10 < L < 50 for the J-Q model), implying that the theory 
correctly captures the mesoscopic physics of competition between the antiferromagnetic and valence- 
bond orders in quantum spin systems. At larger sizes, however, we observe significant deviations 
between the two flows which both demonstrate strong violations of scale invariance. This reliably 
rules out the second-order transition scenario in at least one of the two models and suggests the 
most likely explanation for the nature of the transition in the J-Q model. 



PACS numbers: 64.70.Tg, 05.30.-d, 05.50.+q, 75.10.-b 

The concept of the deconfined critical point (DCP) [T|- 
[5] was developed for understanding quantum transitions 
in two dimensions (2D) between phases characterized by 
different broken symmetries. The key feature of DCP is 
the emergence of fractional degrees of freedom (spinons) 
and gauge fields at the critical point (cf. [1]). Potentially, 
the DCP scenario has a broad range of applications rang- 
ing from quantum phases transitions in lattice models 
and magnets to normal-superfluid transitions in multi- 
component charged superconductors, etc. [IMSl HI IS]- 
Ultra-cold atoms in an optical lattice is another promis- 
ing system where DCP can be tested experimentally [7]. 

A hallmark of the theory is a conjecture that the DCP 
universality class is captured by the 3D classical DCP ac- 
tion involving two complex-valued matter fields, "00=1,2, 
describing spinons coupled to a vector gauge field [T|- 
[31 [5J . Depending on the symmetry group of the underly- 
ing quantum system — global U(l) or global SU(2) — the 
DCP action features the following symmetry in terms of 
its two components: either the Z2 symmetry between two 
spinon fields and the U(l)xU(l) symmetry associated 
with the individual phases of tp a or an enhanced SU(2) 
symmetry between the spinon fields. However, flowgram 
studies of the typical U(l)xU(l) [8 and SU(2) .9] DCP 
actions revealed generic runaway flows consistent with 
weak first-order transitions. 

The initial work focused on microscopic models of 
the superfluid to solid quantum phase transitions first 
claimed the observation of the second-order U(l)xU(l) 
transition [TO] but severe violations of scale invariance 
revealed in the subsequent analysis all but ruled it out 
[TT] , Similarly to the U(l)xU(l) case, early studies of 
the anti-ferromagnetic SU(2)-symmetric J-Q model [TST 
fT4] suggested that the Neel phase transforms into the 
valence bond solid (VBS) in a continuous fashion, while 



subsequent work |151 116j revealed violations of scale in- 
variance. It is important, however, that, up to linear 
system sizes of few hundred sites, the J-Q model clearly 
demonstrates an emergent U(l) symmetry and its run- 
away flow remains rather weak, leaving room for specu- 
lations about the second-order DCP scenario [IB] . 

In this Letter, we perform a direct quantitative com- 
parison of critical flows in the J-Q and the 3D SU(2)- 
symmetric discrete NCCP 1 models. The rationale be- 
hind our study is as follows. Slow run-away flows in 
both models suggest that, independently of the order of 
the transition, the DCP theory in general, and the 3D 
SU(2)-symmetric discrete NCCP 1 model, in particular, 
capture the essence of the quantum phase transition at 
least at intermediate scales of distances. And we indeed 
find that the winding- number flowgrams 8 , 9; of the two 
models can be collapsed in a significantly large region of 
linear system sizes (up to L rj 75 for the J-Q model), 
proving the hypothesis. At larger sizes we observe sig- 
nificant deviations between the two flows which preserve 
their runaway character. The most conservative conclu- 
sion then is that at least one of the two models does not 
feature the second-order criticality, with the straightfor- 
ward interpretation being that both models feature weak 
first-order transitions. 

J-Q and DCP models. The SU(2)-symmetric J-Q model 
describing s = 1/2 spins on a square lattice has been 
analyzed in Ref. [12] : 

h = jJ2"s& -qE - DM - 1) • (i) 

(ij) {ijkl) 

The first sum runs over nearest-neighbor sites (ij) and 
represents the standard Heisenberg model with J > 0. 
The second sum runs over the corners of plaquettes 
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V = (ijkl) such that ij and kl form two parallel ad- 
jacent horizontal or vertical links and adds the four-site 
ring-exchange terms with Q > 0. The model features a 
phase transition between the Neel and an insulating VBS 
phase (its precise nature, however, cannot be determined 
from available system sizes [THE]). While in the Neel 

phase the order parameter S = (S) is linear in the spin 
operator, the VBS long-range order is based on a bilin- 
ear scalar combination of S. Since broken symmetries 
in these phases are different, according to the standard 
Ginzburg-Landau- Wilson paradigm a single phase tran- 
sition between them must be discontinuous. [The actual 
discontinuities, however, should be very weak if one of the 
order parameters is characterized by a significant numer- 
ical smallness far away from the transition, as is the case 
in the VBS phase of the J-Q model where the dimer order 
VD 2 « 1/20 and no signatures of Z4 broken symmetry 
are observed even for largest system sizes [12]]. Thus, if 
a single continuous transition were observed, this would 
be a strong evidence supporting the DCP scenario. 

The DCP is described by the 3D classical two- 
component SU(2) symmetric electrodynamics with the 
emerging U(l) gauge vector-field A [TH3], -^dcp = 
J d 3 x{t\[W - iA]ip\ 2 + ^(V x A) 2 }, where the spinor ip 
consists of two complex fields ip — (ip\, ipz). According to 
the mapping, the Neel vector S = -^n , where n is given 
by 



ft = ip*aip, 



(2) 



with a standing for the Pauli matrices. With the NCCP 1 



fixed-modulus constraint [T], \ipi 



1, one ob- 



tains n 2 — (\ip\\ 2 + |"02| 2 ) 2 = 1 and n + = n x + in y — 
2'/'iV ; 2 implying that the azimuthal angle of n is the rel- 
ative phase of the spinon fields, ip = <p% — ip±, where 
ip a ~ exp(vp a ), a = 1,2. 

The lattice version of the DCP action on a simple cubic 
lattice [H Ej is: 

#dcp = "* E (€A J e^<^+c.c.) + 

(ij), a 
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(3) 



where the gauge field A^ is oriented along the bond (ij) 
from site j to site i, and V x A is the lattice curl-operator 
evaluated on elementary plaquettes V . The effective con- 
stants (t, g) relate in some way to the parameters of the 
J-Q model ([!]). Below we will present evidence that 
g = 1.1 and t = 0.8822(4) provide the closest descrip- 
tion of the J-Q model with J/Q ss 0.04 up to a linear 
size L ~ L* = 75. 

Dual variables. In Ref. [9] , the statistics of the model |3| 
has been reformulated in terms of the dual variables — 
integer bond currents J^ which obey the Kirchhoff con- 
servation laws. Accordingly, the partition function of the 



DCP action Hdcp ^ can be represented as 
Z = J dA Z(Wx,W 2 ) ■ 



exp 



(Stpi 



Wi,W 2 

+ A ]-W 1 +i (6<p2 + A ) ■ W : 
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where Aq stands for the q = harmonic of the gauge field 
defined on the lattice with periodic boundary conditions, 
W a are windings of the bond currents J^ , and S(p a stand 
for the Thouless boundary phase twists of the spinon- 
field phases ip a . By definition, Z(W\, W2) is the partition 
function in a given winding number sector. The integra- 
tion over Aq yields the constraint Wi + W2 — so that 
Z = J2w Z(W,-W) e*v(i8lp-W)wth8lp = 8ip x -5lp 2 . 
The stiffness of the S-vcctor field is found from 
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It is important that at the critical point the scaling be- 
havior of winding numbers is characterized by (W 2 ) — 
0(1) so that ps oc 1/L. In the ordered Neel phase 
(W 2 ) oc L and the stiffness is finite, ps — O(l). 

Our simulations of the J-Q model ([l]) are based on 
the path-integral representation for the partition function 
with periodic boundary conditions in the imaginary time 
< t < /3, where (3 denotes the inverse temperature 
(in both cases we employ the Worm Algorithm approach 
[18] ; simulations of the DCP action were performed as 
described in Ref. [H])- Accordingly, the spin stiffness pjQ 
with respect to the Thouless phase-twist can be expressed 
in terms of the spin world-line windings W' x , W' y along the 
spatial directions x and y , respectively: 



PJQ = Y p [«^) 2 > + ^ W vf)] 



(6) 



In order to compare the two models at the transition 
point we also need to fine-tune the /3/L ratio for each sys- 
tem size L in order to reach the space-time symmetry in 
the J-Q model. We achieve this by defining a space-time 
symmetric winding in the time direction, W' T = S z 

(in the basis where S z — S z — ±1/2 is diagonal), and 
requiring that its mean-square fluctuations coincide with 
{{W' x ) 2 ) = {(W'y) 2 ). We note that W' T is defined without 
the factor of 2 (cf. Eq.(4) of Ref. [IB]). Such definition 
guarantees that fluctuations of W' T proceed in the same 
way as the spatial windings do — in increments of ±1. 

Thus, if parameters of both models (IT]) and ^ are 
kept at the critical point J/Q « 0.04 [T2] and t = t(g) 
(below the bicritical point) [9], the universal values of 
the winding number fluctuations in both models Rjq — 
(l(W x ) 2 + (W') 2 + (W;) 2 ]) ~ and R = ([(W x ) 2 + 

(W v ) 2 + (W z ) ]) ~ O(l) must coincide provided J-Q and 
NCCP 1 models have the same fixed point. 
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FIG. 1: (Color online) Optimal ratio /3(L)/L vs L, with the 
numerical data represented by dots obtained at the pseudo 
critical points defined in the text. Solid red line is the fit 
by D + Aexp(-BL) and the dashed-dotted blue line is the 
fit by D + B/L, with the dashed black line representing the 
asymptote fi/L = D = 0.4270 ± 0.0005 corresponding to the 
space-time symmetry of the J-Q model. 
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FIG. 2: (Color online) Size-dependent transition points 
(J/Q) C (L) of the J-Q model with the j3/L ratios as in 
Fig. [I] The inset shows the pseudo-transition points t c (L,g) 
for g = 1.1 in the DCP model pj. Extrapolation of both 
curves to the L — > oo limit provides estimates of the ther- 
modynamic transition points: (J/Q) c ~ 0.0451 ± 0.0004 and 
t c (g = 1.1) = 0.8822 ± 0.0004. 

Finite size analysis. Simulations of both models have 
been conducted for a sequence of linear sizes using ex- 
actly the same definition of the pseudo-transition point 
in a finite size system, according to the flowgram method 
[H H] . Specifically we tuned model parameters so that 
the ratio of statistical weights of configurations with and 
without windings, J 7 , equals the same constant of order 
unity. We have chosen T = 0.55 because it offers the 
smallest deviations from the space-time symmetry in the 
J-Q model at large L, as shown in Fig. [I] The values of 
the parameters at the pseudo-transition points for both 
models are presented in Fig. [2j 

The universality of scaling behavior is characterized by 
a unique function R = 1Z{J-) in the thermodynamic limit 




FIG. 3: (Color online) Flowgrams of the J-Q (red line) and 
the DCP models (for several values of g). 

L — > oo, [3 L, i.e. for fixed T = 0.55 one expects that 
R(L) curves saturate to the same value even if they de- 
viate from each other at finite L. To see if this is indeed 
the case we have measured Rjq vs L and R vs (L, g) for 
several values of L (from L = 4 to L = 36 for the DCP 
model and from L = 6 to L = 196 for the J-Q model). 
Figure [3] shows the family of DCP flowgrams R(L) for 
several values of the interaction constant g. It also shows 
the flowgram Rjq(L) for the J-Q model. It is immedi- 
ately clear that the values of i?-curves overlap and all by 
itself this is an evidence that DCP theory captures the 
physics of the transition point in the J-Q model. This 
crucial aspect as well as that all the curves feature diver- 
gence with L, in violation of the scale invariance hypoth- 
esis for both models will become more evident below. 

As discussed earlier in Ref. [S], the family of DCP 
flowgrams can be collapsed on a single master curve 
by rescaling system sizes as L — > C(g)L, where C(g) 
is found as a variational distance scale for each value 
of g. This collapse implies that properties of the DCP 
model at coupling strength g — g\ and length scale 
L = L\ are essentially the same as at g — gi and 
L = L 2 = LiC(g\)/C(g2), provided L is larger than some 
microscopic size w 6. Figure [4] shows the quality of the 
data collapse procedure as well as the master curve which 
emerges from it. It also shows the flowgram of the J-Q 
model with rescaled distance L — > CjqL. The value of 
Cjq has been adjusted in order to achieve the best over- 
lap with the DCP-master curve. Note that the freedom of 
choosing Cjq is equivalent to shifting the Rjq curve hor- 
izontally as a whole (in the log L-scale) , i.e. the curve's 
shape remains preserved. It means that the rescaling 
procedure is not supposed to result in the same slope at 
the crossing point between the two flows unless they have 
some common origin. As can be seen, the two curves co- 
incide with each other at length scales 10 < L < 50 (in 
terms of "bare" J-Q model sizes) before they start sig- 
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FIG. 4: (Color online) Flowgrams from Fig.[3]are collapsed by 
rescaling system sizes as L' = C(g)L for the DCP model (this 
amounts to the horizontal shifts of the curves) and L' = 6.8L 
for the J-Q model. Green dotted line shows the master 
curve fit by the A + B(L') a function with A = 0.463, B = 
0.00823, a = 0.437. The lower (orange) dot on the .R'-axis 
indicates the universal value -fi<)(4) ~ 0.475 for the 0(4) uni- 
versality class (g = case). The upper (red) dot on the R' 
axis corresponds to the universal value R' c 
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acterizing the 0(3)-universality. Inset: the rescaling function 
C{g) such that C(0.3) = 1. 
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FIG. 5: (Color online) J-Q flowgram from Fig. [4] is shown 
together with the DCP g = 1.1 flowgram demonstrating the 
best overlap between the two models. The dotted line shows 
the master curve. The vertical arrow indicates the scale (L* = 
75 for the J-Q model ) starting from which the flows diverge 
significantly. The dots on the R' axis mark the 0(4) and 0(3) 
universal values as in Fig. [4j with the dashed horizontal line 
for the 0(3) asymptote. 

nificantly diverging from each other at L > L* = 75. It 
is also important that the J-Q flow starts from the 0(4) 



universal value R' r 



0(4) 



0.475 rather than from the 0(3) 



universality characterized by i?' f 



o(3) 



0.583 as one would 



expect from the classical Heisenberg model, see Fig. [5] 
Finally, as Fig. [5] clearly shows, the J-Q flow runs past 



the 0(3) universality at L > L* . 

Conclusion and discussion. Our key finding is that the 
physics of the transition point between the Neel and in- 
sulating phases in the J-Q model is indeed captured by 
the DCP model up to a large length scale L* = 75. 
At small sizes the flows of R and Rjq start from the 
universal value characterizing the 0(4) universality class 
i^O(4) ~ 0.475. This very fact is a strong indication that 
spinons emerge as dominant degrees of freedom in the 
J-Q model already at length scales L < 8 (in agreement 
with the observed U(l) symmetry of the VBS order pa- 
rameters [12] )■ However, the divergence of the flows at 
L > L* unambiguously excludes the possibility that the 
J-Q model and the DCP action share the same criticality 
in the thermodynamical limit. 

As shown in Ref. [9], the run-away flow of the DCP 
master curve ends up in the first order phase transition 
(detectable at g s» 1.65 for sizes L ~ 30 — 36). [The 
rescaling function C (g) shown in the inset in Fig. [4] is 
a smooth function defined on g > 0. It has no features 
indicating the presence of the tricritical point at some 
g — gtr > 0]. This explains why the J-Q and DCP flows 
ultimately depart. Given the data, there are two possi- 
bilities for the ultimate fate of the J-Q flow: either the 
first order transition or some unknown universality at 
larger values of R' . The fact that both models follow the 
same flow at L < L* and both violate the scale invari- 
ance hypothesis at large length scales strongly favors the 
first possibility — while showing quasi universal behavior 
at intermediate L the two models deviate from this uni- 
versality when the system size is approaching the size of 
the first-order nucleation bubble which does not need to 
be the same in different models. 
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